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In  recent  years,  many  experimental  efforts^"^^  have  been 
made  to  achieve  rapid  plasma  heating  by  the  use  of  intense  relativistic 
electron  beams.  Typically,  a beam  of  50-200  nsec  duration  is  injected 
into  a plasma  or  neutral  gas  column  magnetized  by  a uniform  magnetic 
field  Bjj.  In  such  experiments,  the  most  commonly  used  diagnostic  tool 
to  measure  the  plasma  energy  is  either  a magnetic  probe  or  a diamagnetic 
loop.  To  calculate  the  plasma  perpendicular  energy  density  from  the 
probe  or  loop  measurements,  a sharp  boundary  static  pressure  balance 
model  is  commonly  assumed.  For  example,  if  is  the  measured  (para) 
magnetic  field  variation  between  the  hot  plasma  and  a metallic  wall  of 
radius  r then  from  the  pressure  balance  equation, 

W 

8^,  ♦ 8,tW_  = (8^  + 68^)2  (1) 


where  B^^is  the  magnetic  field  inside  the  hot  plasma,  and  invoking 
conservation  of  magnetic  flux,  can  be  expressed  as 


eq 


z’ 


for  6B,  « B„ 
z 0 


(2) 


where  r^^  is  the  hot  plasma  radius  at  equilibrium. 

However,  static  models  are  not  always  valid.  Under  some 
conditions  (discussed  below),  the  establishment  of  a plasma  equilibrium 
after  rapid  beam  energy  deposition  involves  relatively  slow  mass  motion 
of  the  plasma;  consequently,  on  the  time  scale  of  the  experiment,  the 
equilibrium  state  indicated  by  Equation  1 cannot  be  reached. 


Furthermore,  the  measured  6B^  under  such  conditions  may  be 
the  amplitude  of  a magnetosonic  wave,  instead  of  the  depth  of  a mag- 
netic well.  Therefore,  an  alternative  model  is  needed  when  the  validity 
of  the  static  model  becomes  questionable.  . 
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Physically,  two  limiting  regimes  can  be  distinguished--the 
resistive^^’^^  regime  and  the  reactive^^  regime.  If 

T^.  and/or  T^«  (3) 

where  is  the  ion  collision  time,  x^^^  is  the  magnetic  diffusion 

time,  and  x^  = ''A  Alfven  speed), 

the  system  is  in  the  resistive  regime  and  is  characterized  by  magnetic 
12-15 

diffusion  processes  . The  correlation  between  (the  nonoscillatory) 
and  a given  diamagnetic  current  in  this  regime  has  been  discussed 
by  Guillory  and  Bailey^^,  and  Striffler  and  Kapetanakos^^.  In  the  re- 
active regime, 

1 

""  ■'a  (<)  ' 

< 

the  system  is  characterized  by  gross  plasma  motions,  namely,  magneto-  | 

sonic  oscillations^®’  The  latter  regime  applies  to  plasmas  with 
sufficiently  high  temperatures.  Recently,  oscillatory  magnetic  signals 
I have  been  reported® and  identified®’^ as  magnetosonic  oscilla-  | 

j tions.  With  improved  heating,  future  experiments  are  likely  to  be  in  i 

' this  regime  also.  Hereafter,  we  will  be  concerned  with  this  regime 

i 

only. 

It  is  important  to  distinguish  two  cases  according  to  experi- 
mental conditions.  Case  1:  the  hot  plasma  is  surrounded  by  vacuum 
or  by  a neutral  gas;  Case  2:  it  is  surrounded  by  a cold  plasma.  For 
case  1,  B^  between  the  plasma  and  the  wall,  though  oscillatory,  is 
spatially  uniform  (Fig.  la),  thus  the  static  model  gives  a reasonable 
estimate  of  provided  one  takes  the  time  averaged  5B^  and  determines 
J Tgq  consistently.  For  case  2,  on  the  other  hand,  the  probe  or  loop 

(immersed  in  the  cold  plasma)  will  measure  a magnetosonic  wave  (Fig.  lb) 
generated  by  the  beam  energy  deposition  and  the  correlation  between  6B 
and  is  determined  by  the  dynamics  of  the  plasma,  rather  than  by 
the  condition  of  static  pressure  balance.  The  implication  of  6B^  in 
this  case  is  thus  qualitatively  different  from  that  of  case  1.  For 
example,  roughly  speaking,  the  wall  would  have  no  effect  on  the  peak 
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probe  signal  if  i-fp  and  it  has  no  effect  on  the  peak  loop 
signal  if  r^<  j Tp,  where  rp  is  the  plasma  radius,  r^  is  the  probe 

position  for  loop  radius  and  we  have  assumed  constant  wave  speed. 

Clearly,  if  the  above  conditions  are  satisfied,  W given  by  Eq.  (2), 

2 

which  scales  as  r^,  is  Inapplicable. 

In  most  experiments, the  plasma  is  preformed  by  a 
discharge;  it  is  reasonable  that  in  these  experiments  case  2 prevails. 

To  the  authors'  knowledge,  a proper  theoretical  model  to  interpret  probe 
or  loop  data  for  case  2 in  the  reactive  regime  has  hitherto  been  lacking, 
although  there  are  ample  experimental  evidences  that  this  regime  has 
indeed  been  observed.  In  the  following,  we  develop  a dynamical  model 
applicable  to  this  regime. 


Formulation 

We  make  the  following  assumptions:  (1)  3/30  = 3/3z  = 0; 

(2)  * Vg2  “ Oi  (3)  initially  uniform  plasma  (0=0^^)  maintaining 

quasineutrality;  this,  together  with  the  previous  two  assumptions, 

implies  v^^  = v^^  = v^.;  (4)  v^  « c;  (5)  neglecting  electron  inertia 

and  ion  pressure;  (6)  neglect  electron  and  ion  collisions;  (7)  neglect 

2 -1 

radial  thermal  conduction— this  requires  "^A* 

Pg  is  the  electron  Larmor  radius  and  Vg  is  the  electron  collision  fre- 
quency; (8)  neglect  wall  effect;  and  (9)  isotropic  electron  velocity 
distribution.  The  system  is  then  described  by 


(—  * V ■ fp 

3t  3r  ^ 


(5) 


^r  2 — ^e 

^ ^ 3r  ® 


(6) 
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(7) 


b 

r 3r  ^ 3t  ^ 

(8) 

i_  17  + -L  1^7  7 7j  - 0 

3t  r 3r  ” 

_5  _5 

(9) 

V “ 1" 

3t  3r  ® ^ 

(10) 

where  overbars  represent  normalized  quantities  defined  as  follows: 
t = t/T^,  ? i r/rjj,  rT  = n/n^,  Fg  = Sirnklg/Bp,  = B^/B^, 

"^0  - \ - V/''A’  ^b  = ^b  ^^V^o 

Qjj  is  the  rate  of  beam  energy  deposition  per  unit  volume.  We  let 


Ob' 

where  normalized  beam  duration,  the  free  parameter 

1 . 

0£  * - / 2rrrdr  / dt  Qj^,  is  the  space  and  time  integrated  beam 

■'o  •'O 

energy  deposition  scaled  to  the  total  external  magnetic  field  energy 
within  the  beam  volume  (r  £ r^^)  and  s is  a steepness  parameter  to 
account  for  the  spatial  inhomogeneity  of  beam  energy  deposition. 

Note  that  the  only  parameters  to  specify  are  those  contained  in 
Qjj  (i.e.,  s,  Tjj,  and  8^).  All  other  parameters  such  as  B^,  n^,  r^^, 
etc.,  are  scaled  out  of  Eqs.  (5)-fl0)  through  normalization. 


Lx.  ^ sin  , if  r < 1,  t < 1 


0 , otherwise 
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Solutions 


The  radial  one-dimensional  evolution  of  the  system  is  computed 

from  Eqs.  (5) -(10)  using  a hydromagnetic  code  which  employs  a two-step 

Lax-Wendroff  flux  corrected  transport  (FCT)  algorithm. 

» 

Results  of  our  calculations  are  schematically  presented. 

Figure  2a  shows  typical  radial  profiles  of  Figure  2b  plots 
typical  probe  signals  [58^(7^)]  and  loop  signals 

<S?(r)  H dP  2ttP  58  (P) 

® J Q * 

as  would  be  measured  by  a probe  or  loop  with  7^  =•  1.2.  If  a conducting 
wall  is  present,  wave  bouncing  will  predictably  produce  a second, 
third  peak,  and  so  on. 

Figures  3a  and  b plot  the  peak  probe  signals  (58^)  and  loop 
signals  (6$^)  versus  the  total  energy  deposited  (S^).  Each  family  of 
curves  is  generated  by  the  parameter  We  note  that  for  the  same 
amount  of  energy  deposited,  a faster  rate  of  deposition  (i.e.,  smaller 
T^)  results  in  larger  wave  amplitude  (see  the  68^  curves),  and  stronger 
dispersion  or  nonlinearity  (see  the  bending  of  the  5^^  curves).  Such 
differences  can  be  distinguished  in  the  present  dynamical  model,  but 
not  in  the  static  model.  Since  58^  is  a local  quantity  and  5i*^  is  a 
global  quantity,  7^^  has  less  effect  on  the  latter. 

All  calculations  presented  so  far  are  for  the  steepness 
parameter  s=2.  We  have  varied  s in  our  runs  and  found  that  the  s*l 
case  gives  slightly  (<15%)  higher  and  the  s=3  case  gives  slightly 
(<10S)  lower  68^.  The  effect  of  s on  5$^  is  much  weaker. 

Because  of  the  normalized  parameter  system  employed  here, 

the  data  shown  in  Figs.  3a  and  b have  general  applicability.  As  an 

example,  consider  the  Physics  International  experiment  recently  reported 

by  Prono,  et  al.^®  Their  high  v/y  beam  heating  scheme  produced  an 

L9  3 

impressive  plasma  energy  density  as  high  as  10  eV/cm  . There  is 
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little  question  that  the  experiment  was  in  the  reactive  regime.  The 

static  model  that  they  used  to  interpret  the  probe  signal  is  appropriate 

for  the  situation  depicted  in  Fig.  la.  If  the  situation  was  more  like 
21 

Fig.  lb  instead,  then  the  dynamical  model  should  have  been  used.  In 
any  case,  it  is  interesting  to  compare  their  interpretations  with 
those  given  by  the  dynamical  model.  Using  their  published  data,  we 

obtain  7 ^ 1.2,  Tu  1-2,  and  68^  0*1»  thus  from  Fig.  4a, 

^ ^ ^ 20 
Bp  which  corresponds  to  an  average  plasma  energy  density 

^ 18  3 

of  4-6  X 10  eV/cm  , In  agreement  with  their  interpretations.  We 
point  out,  however,  that  the  agreement  between  the  two  interpretations 
is  coincidental  in  view  of  the  qualitative  differences  between  the 
models  used. 

Finally,  if  either  collisions  (magnetic  diffusion  effect)  or 
thermal  conduction  (cooling  effect)  were  included  in  our  model,  the 
calculated  signals  (68^  and  5^)  would  be  weaker.  Thus,  8^  inferred 
from  the  present  model  yields  a conservative  estimate  of  beam  energy 
deposition. 


The  authors  would  like  to  thank  Drs.  D.  8ook,  D.  A.  Hammer, 
and  C.  C.  Wei  for  helpful  discussions. 
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The  present  fluid  model  only  gives  isotropic  energy  density. 

It  cannot  distinguish  the  Maxwellian  distribution  from,  for 
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Experimentally,  this  question  can  be  resolved  by  measuring  6^  at 
two  radial  positions  outside  the  beam  channel  and  observing  t^eir 
phase  relationship. 
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Director 

Ballistic  Missile  Defense  Program  Officer 
Commonwealth  Bldg., 

Arlington,  Va.  22209 

ATTN:  Dep  BMDPM  Mr.  Julian  Davidson 

Mr.  Archie  Gold 


Director 

US  Army  Ballistic  Research  Laboratories 
Aberdeen,  Maryland 

ATTN:  Mr.  Frank  Niles 

Dr.  J.  Martlneau 

Centre  de  I'Energie 

University  of  Quebec,  INRS  CP  1020 

Varennes,  Quebec,  CANADA 

Dr.  R.  Papoular 

Conm.  Energle  Atomlque  - France 

B,  P,  No,  6 - Fontenay  - Aux  - Roses  - 92 

France 

Dr.  Ralph  Rudder 

Air  Force  Weapons  Laboratory 

LRE  - Kirkland  AFB 

New  Mexico  87117 
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Dr,  David  Smith 

United  Aircraft  Research  Laboratories 

Loo  East  Main  Street 

East  Hartford,  Conn.  06IO8 

Dr,  R.  G.  Tomlinson 

United  Aircraft  Research  Laboratories 
Silver  Lane 

E.  Hartford,  Conn.  06IO8 

Dr.  Robert  Turner 
Applied  Physics  Laboratory 
Johns  Hopkins 
8621  Georgia  Avenue 
Silver  Spring,  Md.  209IO 

Or.  George  Vlaska 
Aerospace  Resea:  ch  Laboratory 
University  of  Washington 
Seattle,  Washington  98IO5 

Dr.  H.  J.  Kunze 

Abteilung  fur  Physik  und  Astronomie,  Ruhr  - Universitat 
L63  Bochum 
Postfach  21L8 
West  Germany 

Dr,  A,  W.  DeSilva/Dr.  R.  C.  Davidson/  Dr.  J.U.  Guillory 

University  of  Maryland 

Department  of  Physics  and  Astronomy 

College  Park,  Maryland  20Jk2 

Dr.  A.  J.  Alcock 
National  Research  Council 
Physics  Division 
100  Sussex  Drive 
Ottawa,  Ontario,  Canada 

Dr.  Martin  C.  Richardson 
National  Research  Council 
Physics  Division 
100  Sussex  Drive 
Ottawa  2,  Ontario,  CANADA 

Dr.  Keith  Boyer 

Los  Alamos  Scientific  Laboratory 
P.O.  Box  1663 

Los  Alamos,  New  Mexico  875^ 

Dr,  John  W.  Delber 
Dr.  Ron  Rehm 

Cornell  Aeronautical  Laboratory 
LL55  Genesee  Street 
Buffalo,  New  York  1L221 
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Dr.  D,  F.  Dubois 
Hughes  Research  Laboratories 
3011  Malibu  Canyon  Road 
Malibu,  California  90265 

Dr.  A.  G.  Engelhardt 
Hydro-Quebec  Institute  of  Research 
1800  Montee  Ste- Julie,  Varennes 
P.  Q.,  Canada 

Dr.  T.  V.  George 

Westinghouse  Research  Labs.  Pgh. 
Pittsburgh,  Pa.  15235 

Dr.  Martin  V.  Goldman 
Department  of  Astro-Geophysics 
University  of  Colorado 
Boulder,  Colorado  80302 

Dr.  Arthur  H.  Guenther 
Air  Force  Weapons  Laboratory  (S.Y.) 
Kirtland  Air  Force  Base 
Albuquerque,  New  Mexico  87117 

Mr.  Roger  Case 
Air  Force  Weapons  Laboratory 
Kirtland  Air  Force  Base 
Albuquerque,  New  Mexico  87II7 

Dr.  Eric  D.  Jones 
Sandia  Laboratories 
Div.  5214 
P.O.  Box  5800 

Albuquerque,  New  Mexico  87II5 

Dr.  Moshe  J.  Lubin 
University  of  Rochester 
Laboratory  for  Laser  Energeties 
Mech.  and  Aerospace  Sciences,  Dept. 
Rochester,  N.  Y.  14627 


Dr.  Philip  Mallozzi 
Battel le  Memorial  Institute 
Dept,  of  Physics 
505  King  Avenue 
Columbus,  Ohio  43201 

Dr.  James  Shearer 
Lawrence  Livermore  Lab. 

P.  0.  Box  808 

Livermore,  California  9^551 
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Dr.  Donald  W.  Kerst/  Dr.  R.  S.  Post 
University  of  Wisconsin 
Sterling  Hall 
Department  of  Physics 
Madison,  Wisconsin  53706 

Prof,  J.  G.  Hirschberg,  Jr, 

University  of  Miami 
Department  of  Physics 
Coral  Gables,  Florida  33124 

Dr.  L.  P.  Bradley 

Lawrence  Livermore  Laboratory 

P.O.  Box  808 

Livermore,  Calif.  9^551 

Prof.  Frank  Chen  / Prof.  M.C.  Luhmann 
Department  of  Engineering 
University  of  California 
Los  Angeles,  Calif.  90024 

Prof.  Alfred  Wong  / Prof,  J.  Dawson 
Department  of  Physics 
University  of  California 
Los  Angeles,  Calif.  90024 

Dr.  Stephen  0.  Dean 
Division  of  Research 

Energy  Research  and  Development  Administration 
Washington,  D.C.  20545 

Dr,  Roy  Gould 
Cal.  Inst.  Tech. 

Pasadena,  Calif,  9HO9 

Dr,  J.  Decker/  Dr,  W.  F.  Dove 
DCTR 

Energy  Research  and  Development  Administration 
Washington,  D.  C.  20545 

Mr.  W,  C.  Gough  / F.  R,  Scott 
Electric  Power  Research  Institute 
Washington,  D.C.  20545 

Dr.  S.  J.  Buchsbaum 
Bell  Telephone  Labs 
P.  0.  Box  262 
Murray  Hill,  N.  J.  07904 

Dr.  D.  J.  Rose 

Oak  Ridge  National  Lab. 

P.  0.  Box  X 

Oak  Ridge,  Tenn.  37830 
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Or.  A.  E.  Ruark 
7952  Orchid  St.  N.  W. 
Washington,  D.C.  20012 

Maxwell  Laboratories 
9244  Balboa  Avenue 
San  Diego,  California  92123 
ATTN:  Dr.  A.  C.  Kolb 

Dr.  Peter  Korn 
Dr.  V.  Fargo 
Mr.  J.  Shannon 

Dr.  A.  F.  Haught 
39  *ox  Den  Road 
Glastonburg,  Conn.  06033 

Dr.  F.  J.  Fader 
197  Foote  Road,  S. 
Glastonburg,  Conn.  060T5 

Dr.  Fred  Schwirzke 
Naval  Postgraduate  School 
Monterey,  Calif.  92940 

Dr,  M.  B,  Gottlieb 
Forrestal  Research  Center 
P.  0.  Box  451 
Princeton,  N.  J.  08540 

Dr.  F.  L.  Ribe 
P,  0.  Box  1665 

Los  Alamos,  New  Mexico  87544 

Dr,  T.  K.  Fowler 
Lawrence  Livermore  Lab. 

P.  0.  Box  808 
Livermore,  Calif.  9455I 

Dr.  Herman  Postma 
Bldg.  9201-2  336 

T-12 

Oak  Ridge,  Tenn.  75830 

Dr,  J.  Rand  McNally 

105  Norman  Lane 

Oak  Ridge,  Tenn,  73830 

Dr,  Wulf  Kunkel 
Dept,  of  Physics 
University  of  California 
Berkeley,  Calif,  94720 

Dr,  Harry  Dreicer 
P.O.  Box  1653 

Los  Alamos,  New  Mexico  87544 
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Dr.  Franz  Johoda 
P.  0.  Box  1665 

Los  Alamos,  New  Mexico  875^^ 

Dr.  Ralph  Lovberg 
14.714-14.  Panorama  Drive 
San  Diego,  Calif.  92116 

Dr.  Herbert  W.  Friedman 
Avco  Everett  Res.  Labs. 

2585  Revere  Beach  Parkway 
Everett,  Mass.  02lk<^ 

Cornell  University 
Ithaca,  New  York  11^850 
ATTN:  Dr.  David  Morse 

Dr.  R.  N.  Sudan 
Dr.  C.  B.  Wharton 
Dr.  Hans  Flelschman 
Dr.  J.  A.  Nation 

Dr.  Joseph  Fergusen 
Dept,  of  Physics 
Mississippi  State  University 
State  College,  Miss.  39762 

Dr.  Robert  J.  Mackln,  Jr. 

2622  N.  Holliston 
Altadena,  Calif.  9^001 

Dr.  David  R.  Bach 

8 Heatherldge 

Ann  Arbor,  Mich.  k-BlCh 

Dr.  Richard  Hall 
1851  111  E.  N.  E. 

Bellevue,  Washington  980024- 

Dr.  Lawrence  Lldsky 
17824-  Washington  St. 

Newton,  Mass.  02166 

H.  Ahlstrom  / L.  Steinberg 

Dept,  of  Aeronautics  and  Astronautics 

Univ.  of  Washington 

Seattle,  Wash.  98IO5 

Dr.  David  Koopman 
Inst.  Fluid  Dynamics 
University  of  Maryland 
College  Park,  Md.  207*4-0 

Columbia  University 
New  York,  N.  Y.  10023 
ATTN:  Dr.  Robert  Gross 

Dr.  S.  P.  Schleslnger 
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